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ABSTRACT 



The minimal supersymmetric standard model (MSSM) with universal soft breaking terms 
has a negligible effect on the rare c — > wy decay rate. We therefore study a general model 
with non-universal soft breaking terms constrained by vacuum stability and the experimental 
upper bound on the mass difference of the neutral charm mesons. It turns out that gluino 
exchange can enhance the standard model rate by up to two orders of magnitudes; in contrast, 
the contributions from charginos and neutralinos remain at least an order of magnitude below 



the QCD corrected standard model rate. 



1. Introduction 

Flavour changing neutral currents (FCNC) occur in the standard model only at the 
loop level where they are suppressed by the GIM mechanism. The corresponding processes 
are rare and suitable as probes for new physics, in particular for the popular low energy 
supersymmetry (SUSY) which has various new sources of flavour violation. Studying the 
FCNC processes may reveal the flavour structure of the soft supersymmetry breaking terms 
and shed light on the mechanism of SUSY breaking. There is a renewed interest in this 
issue, because models with effective supergravities U generally do not lead to universal soft 
breaking terms and tend to have large FCNC rates. 

In this note we explore the sensitivity of the c — > wy decay rate in the framework of the 
minimal supersymmetric standard model (MSSM) with unbroken R-parity and arbitrary 
non- universal soft breaking terms. This process was considered before in ref. 0; however, 
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only gluino exchanges were included and the constraints from vacuum stability were not 
applied. In addition, we take into account the QCD corrections within the standard model, 
which are seizable, and discuss the QCD corrections in the MSSM. Motivated by recent 
proposals for probing c — > wy in specific hadronic observables || |], ||, which are largely 
free of the 'uninteresting' long distance contributions, we analyze this decay in more detail. 

In the standard model, the c — > u-y amplitude is proportional to J2 q =d s b KgKj9 m <j a ^ the 
one-loop level and the corresponding branching ratio ~ 1CT 17 is negligibly small as the GIM 
mechanism is extremely effective in charm decays (the intermediate d and s quarks are very 
light compared to the W). The most important QCD corrections arise from the diagrams 
shown in Fig. 1. They drastically lift the GIM suppression; instead of power-like it becomes 
logarithmic [|7|. The effective local interaction can be written as @ 

£e ff = - 4 -^v; s v us [or o 7 + c> 7 O'] (1) 

7 = j^m c r«vi(l + 75 )c , 0' 7 = j^m c ua^(l - 75 )c F^ 
(c 7 )ff ~ ^A C2 ( mc ){f[ml/ml] - f[m 2 d /ml]} ~ -(0.007 + 0.020 i)[l ± 0.2] 
(c' 7 ) 5M ~0 

for m c = 1.5 GeV and a s (m c ) = 0.31, with functions / and C2 given in 0. The corresponding 
branching ratio in the standard model 
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is still small and potentially sensitive to various extensions beyond the standard model. 
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Figure 1: The most important QCD corrections to c — > w~) are induced by the four-quark 
operator (9 2 at the two loop level. Only one of the representative diagrams in the standard 
model is shown. 

The main generic problem for probing the c — > wy transition in hadron decays is that 
the most straightforward observables such as the decay rate D — >• py are largely dominated 
by long distance contributions and are independent of new physics ||. However, there are 
observables which are dominated by the interesting short distance physics. For instance, 
the disturbing long distance contributions largely cancel in the difference of the decay rates 
Br(D° — > p°7) — Br(D° — > uy) as shown in ||. Other possibilities to look for c — > uy are 
the B c — > Z?*7 decay Q and the violation of the equalities Br(D+ — ► K* + >y)/Br(D+ 
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p+ 7 ) = tan0 c |], Br(D° -> p°^)/Br{D° -> iP° 7 ) = tan# c /2 |, |. Similarly, the isospin 
violating quantity r(.D + — > p + 7)/2r(D° — >■ p°7) — 1 || might help to disentangle the dif- 
ferent contributions. These observables are suitable for testing new scenarios that could 
significantly enhance the c — ► wy rate, and are of particular interest in the context of SUSY 
flavour problem. 




Figure 2: One-loop digrams for c — > wy decay in MSSM. Here qi=i r . } e, Xj=i,2 an d Xi=x..4 
denote mass eigenstates for squarks, charginos and neutralinos. The photon is attached to 
any charged line. 



2. The MSSM amplitudes 

The one-loop diagrams for c — > wy transition in MSSM are shown in Fig. 2. The <7i=i,..6, 
Xj=i 2 an d X°=i .. 4 denote mass eigenstates of squarks, charginos and neutralinos, respectively. 
We calculate the Willson coefficients at the high energy {my/) scale for different contributions 
using the Feynman rules and notation presented in |10| 
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with e d = —1/3, e u = 2/3 and 
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The expressions for c' 7 are analogous, but they are multiplied by an overall factor m u /m c 
and m c is replaced by m u in @ . The matrices Z ± are the mixing matrices of chargino mass 
eigenstates Xj=i,2 with weak eigenstates (W + ,H + ) and the matrix Z^ describes the mixing 
of the neutralino mass eigenstates x?=i .. 4 with weak eigenstates (B, A 3 , H®, H®) [F0| - The 



functions -Fi,..,4 are given in Appendix B of ||TT|. The matrices Zu t u diagonalize the squark 



mass matrices and depend on a large number of unknown soft breaking parameters in the 
general MSSM. 

The number of these parameters entering the c — > wy amplitude is significantly reduced 
in the mass insertion approximation, where the expressions (S) are expanded to the first 
order in the physical squark mass splitting around the common average squark mass rriq 
pj . i.e. 
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This approximation works well |13] unless there are strong cancellations or large flavour 
violations (for instance a light sbottom with a mass around 4 GeV [Q could induce large 
flavor- violating effects. However, its contribution to c — > wy is negligible Q). 

The summation over the intermediate squarks can be expressed in terms of the squark 
mass matrices M 2 ^ and M. 2 ^ in the super-CKM basis |T2|]. In this basis for the squarks, 
the flavour structure of quark-squark-gaugino vertex is the same as the quark-quark-gauge 
boson vertex. We have 



ziizff* = J2 zpz£ = s IJ , Yl z u z u*< 
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In evaluating (Q), we neglect the masses of the quarks except m c because the magnetic 
transitions O-j and 0' 7 (|1]) are of this orderQ. The top quark mass m t does not enter the 



1 The sbottom-chargino exchange gives (c?y 
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<x V* b V ub Z%*Z% F(ml + /ml) m 2 w /m? cc V; b V ub Z}j*Z% 



in the limit m~/m~ + <C 1, where the unitarity of Zu has been applied. 



We neglect the terms proportional to m u /m c in c' 7 . 
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amplitude and m b comes only in combination with a small factor V* b V ub . This approximation 
is supported by the fact that c\ v+ oc V* q V uq m 2 q (||) is about four orders of magnitude smaller 
than (cj)^ 1 (H). Using the relations and the explicit expressions for the squark mass 
matrices^ we getfj 
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The relevant diagrams for the various contributions to C7 (0) in the mass insertion ap- 
proximation are shown in Fig. 3. The expression V(M^)llV occurring in the chargino 
contributions (^) is equal to (M^) LL at the unification scale |TJ|. At the weak scale, this 
relation is valid up to the corrections proportional to quark masses, which can be neglected 
within our approximation. The values of c' 7 and the left-right induced C7 are equal provided 



|(M?))j L |, which is indeed the case for the upper bounds given below 



QCD corrections will not drastically change our main result. As we will see, the contri- 
butions from charginos, neutralinos and left-left induced gluinos are generally smaller than 
those in the standard model. Thus the main effect found in ref. J!]], namely that the cor- 
rections arising from O2 are dominant, remains true for them. We need to include the QCD 
corrections only for the left-right induced gluino contribution because it might give a larger 



sj. 



3 A similar procedure is used in the Appendix of |12| for b 

4 We omitt the terms proportional to (M^,) RR in c' 7) since the best enhancement in general MSSM is 
expected when the upper bound on Amp is saturated only by diagrams with with (M?)^ insertions (see 
Table 1), while (Mf) R 2 R is set to zero. 
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branching ratio than the ~ 3 x 1CT 8 of the standard model. Using Eqs. (37) and (45) of 
ref. [13 1, one finds that the anomalous dimensions of the relevant operators induce a 10% 
decrease of the amplitude ([!]), which is evaluated using ot s = a s {m w ) = 0.12. 



3. The constraints on the mass insertions 

The Wilson coefficients (|7|) depend on the mass insertions (M?)P L , (M~ ) 12 R and (M~ ) RL 
in the super-CKM basis which incorporate the model dependence arising from the mechanism 
of supersymmetry breaking. In MSSM with universal soft breaking terms (i.e. models with 



gauge mediation |15fl) the squark mass matrices are flavour diagonal at the unification scale. 
Flavour mixing is due to the renormalization group evolution down to the weak scale, because 
the Yukawa matrices for up and down squarks are not simultaneously diagonalizable. The 
off-diagonal elements in squark mass matrices at the weak scale are necessarily proportional 



to the product of Yukawa couplings, K iv±^ JLL , 
effect of this constrained MSSM scenario on c - 



(Mp 12 



and the 



«7 decay rate is negligible. 

Supergravity scenarios Jl| generally lead to the non-universal soft-breaking terms and 
FCNC can be used to probe their textures. The mass insertions are taken as free pa- 
rameters in this case, with bounds coming from available experimental data and the con- 
sistency of the model. The mass insertions (M~) 12 can be constrained by saturating the 



upper experimental upper limit on Amp, as done in [ l7| . Taking the recent CLEO result 
(Amp cos5 + ATd sin 5/2)/ (2T 2 D ) < 0.04% |18| and assuming the relative strong phase 6 
between D° — > K + n~ and D° — > K~n + to be small, the upper bounds on mass insertions 
are given in Table 1. 





\(M 2 )l\\/ml 


\{M 2 f 2 R \/m 2 
l(M?)^|/m? 


0.3 
1.0 
4.0 


0.03 m g -/500 GeV 
0.06 m 9 -/500 GeV 
0.14 mq/500 GeV 


0.04 rn^/500 GeV 
0.02 m^/500 GeV 
0.02 m^/500 GeV 



Table 1: Limits on (M?)^ 2 L , {M^) 12 R and (M?.))^, obtained by saturating the upper bound 
Am D < 4.5 x 10~ 14 GeV of CLEO [111 by gluino exchange O . 



Stricter upper bounds on (M?)^, and (M?.)^ 2 L are obtained by requiring the minima of 
the scalar potential in the MSSM not to break electric charge or color and that they are 
bounded from below [|l9j. In particular 



|(M?)i 2 H | |(M?)^ 2 L | r m c 500 GeV 

% , u ' <V3 — ~ 0.0023 . 8 

777- 771- 777,^ 777,^ 

The factor m c in (^) cancels the factor m c in {Ml) l2 R /m c (^), which in turn comes about 
since the chirality flip in the squark state does not have to be accompanied by the chirality 
flip in the external c quark state. 
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Another upper bound 



12 



(M?)i 2 L < max 



(9) 



is due to the relation (M^) LL = ViM 2 ^) LL V^ . Although the constraint on (M|)^ from 
Arrix |17| i s about 30% stronger than the constraint on (M?)J-^, given in Table 1, the upper 
bound (|]) is ineffective due to the large ratios V 12 /V 32 and V 12 /V 13 . 

We have explicitly verified that the upper bounds on {M^) 12 L (Table 1) and (M^) 12 R 
(§) do not saturate the experimental value (Am^) exp ~ 3.5 x 1CT 15 GeV through chargino 
exchange^]. 

4. The results 

We consider low-energy supersymmetry with masses of superpartners in the TeV range. 
We take 250 GeV< m^rrig < 1000 GeV with lower bounds given by the direct searches for 



these states [21]. The masses of charginos, neutralinos and the corresponding mixing matrices 
Z ± , Zn can be expressed in terms of tg/3, /i, mi and m 2 |10], |22|. We take tg(3 and \i in the 
range 2.5 < tgf3 < 30 and 100 GeV < < 300 GeV. We assume mi ~ m 2 and determine 
this parameter by setting the mass of lighter chargino equal to the experimental lower bound 
m x + > 90 GeV (we vary this parameter in the range 90 — 200 GeV). Points of parameter 
space, where the neutralino masses violate bounds from the direct experimental searches, 
are omitted. The technical problem of negative eigenvalues for chargino and neutralino mass 



matrices is solved following the discussion in Appendix A. 3 of [22|. The value a s = a s (mw) = 
0.12 is used. 

The predicted c — > wy branching ratios are given in Fig. 4 in a form of a scatter plot, 
with input parameters varying in ranges discussed above. The mass insertions (M?)^, and 
{Mq)^ rrl are taken at the maximal values of the allowed ranges, given in Table 1 and Eq. 
(|), respectively. The QCD-corrected standard model prediction (JJ, 0) is represented by 
the dot on the left. The remaining dots represent branching ratios arising solely from the 
genuine supersymmetric contributions; as discussed above, QCD corrections hardly affect 
these (0) and are neglected. The rates that would result from the gluino, the chargino 
and the neutralino contributions (|7|) are shown separately. If {M^) 12 R is close to the upper 
value given by the vacuum stability argument (§), the gluino exchange is dominated by 
the left-right mass insertion and can enhance the standard model rate (Q) by two orders of 
magnitudes 

Br[c -> u~f) LR ~ 6 x 10~ 6 (10) 

for rrig ~ rrig ~ 250 GeV. The process c — > w~y is therefore particularly suitable for probing the 
universality of the trilinear soft-breaking A terms, that are responsible for the left-right mass 
insertions. If the dominant supersymmetric contribution comes from the (M-.)]^ insertion 
(when (M?) 12 R is negligible), then the gluino and chargino contributions are comparable 

5 The upper bound (M?)^ (||) gives a Am K which is several orders of magnitudes below (Am K ) exp . 
The upper bounds (M?.)^ 2 L (Table 1) can render values close to, but not above (Am^-) e:cp , for low masses 
of the charginos. The expressions for the chargino contributions to Am are presented in EQ]. 
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in size and in fact add up. In this case, the highest rate of genuine SUSY contribution is 
achieved at small tg/3 and p < 

Br{c -> m 7 ) ll ~ KT 9 (11) 
and would be screened by the standard model contribution (|2|). 



5. Conclusions 

We have examined the sensitivity of the c — - > wy decay rate in the context of the minimal 
super symmetric standard model. Not surprisingly, schemes with universal soft breaking 
terms have a negligible effect on this process. If the universality condition is relaxed, the 
chargino and neutralino contributions remain small and are at least one order of magnitude 
below the standard model rate. However, the gluino exchange can enhance the standard 
model rate of Br(c — > wy) ~3x 1CT 8 by up to a factor of 200, given the various constraints 
on the relevant left-right mass insertions. This is quite a surprising result and makes this 
decay very attractive for flavor studies in supersymmetric models. To understand this drastic 
effect, we recall that the standard model branching ratio is essentially due to a two- loop 
diagrams in Fig. 1, because the one loop-contribution is heavily GIM suppressed. As a 
consequence, one may expect that the GIM non-suppressed two-loop amplitude is enhanced 
by a factor a s {mw) / \otwV* s V us a s {m c )) in a general MSSM, in qualitative agreement with 
our detailed numerical result. 

Although unwanted long distance effects strongly affect the hadronic rates, it is possible 
to find observables where they largely cancel. An important observable is R = [Br(D° — > 
wy)—Br(D° — > p 0r y)]/ Br(D° — > u>j) which can go up to 0(1) in the general MSSM discussed 
here, as compared to a standard model prediction of 6± 15% ||. Although the long distance 
effects are not so well under control ||25|| , a value of R above 20% would be a clear signal 



for new physics. Present experimental upper bounds on the rare decays D° — > p°7 and 



D — > p to are at the level of 10 [G3, 21], while theoreticaly predicted are at a level of 10 
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H. It is expected that BaBar, Belle, BTeV and possibly a tau-charm factory can observe 
these interesting decays. 

Another observable which is relatively free of long distance effects is Br(B c — > B^y). 
While the standard model expectation is 1 x 10~ 8 0, it could reach 1 x 10~ 6 in the MSSM. 
Such an enhancement can be observed at LHC, which is expected to produce 2 x 10 s B c 
mesons with pt(B c ) > 20 GeV. 
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(M?)i 2 L (M?) 




(a) Dominant gluino exchange diagrams. 




(b) Dominant chargino exchange diagrams with flavour index a, b — 1, 2, 3. 




B,A 3 

(c) Dominant neutralino exchange diagrams. 



Figure 3: The one-loop diagrams for c — > try decay in MSSM within the mass insertion 
approximation. Only the diagrams, that give raise to the local operator 7 oc ul&^cr (fj), 
are shown. Quark masses are neglected and only the terms of the order of 0(m c ) are kept. 
The stars denote mass insertions for quarks and squarks. The squark mass insertions M?. ^ 
and primed squark fields q' are represented in super-CKM basis. The dots denote mixing 
of different weak-isospin chargino or neutralino eigenstates. The photon is attached to any 
charged particle. 
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Figure 4: The predicted c — > wy branching ratios for various contributions in MSSM. The 
Q CD-corrected standard model rate ([[], |2|) [0] is represented by the dot on the left. The 
remaining dots represent the rates arising solely from various supersymmetric contributions 
and do not include the QCD corrections (|7|): the gluino exchange (the contributions arising 
from (M?)^ RL and (M?.)^ 2 L are shown separately), the chargino exchange and the neutralino 
exchange (the contributions arising from (Mf )^ RRL and (M?)^ are shown separately). 



The mass insertions {M^) f L and {]^V)lrrl are taken at the maximal values of the allowed 
ranges, given in Table 1 and Eq. (§), respectively. Other parameters are varying in the 
ranges, discussed in the text. 
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